Can the scale factor be rippled? 
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We address an issue: would the cosmological scale factor be a locally oscillating quantity? This 
problem is examined in the framework of two classical 1+1-dimensional models: the first one is 
a string against a curved background, and the second one is an inhomogeneous Bianchi I model. 
For the string model, it is shown that there exist the gauge and the initial condition providing 
an oscillation of scale factor against a slowly evolving background, which is not affected by such 
an oscillation "at the mean". For the inhomogeneous Bianchi I model with the conformal time 
gauge, an initially homogeneous scale factor can become inhomogeneous and undergo the nonlinear 
oscillations. As is shown these nonlinear oscillations can be treated as a nonlinear gauge wave. 

PACS numbers: ll.25.-w, 04.20.Ex, 98.80.-k 



I. INTRODUCTION 

Cosmological constant problem is a subject of numer- 
ous investigations The fact of the matter is that the 
enormous energy of zero point quantum fluctuations of 
matter fields as well as electromagnetic and gravitational 
waves have to produce an enormous universe expansion, 
which is not seen in reality. At the same time, a study 
of conformal fluctuations also has a long history [a-flll] 
and is still far from the completion owing to lack of a 
developed theory of quantum gravity. 

A possible relation between these two problems has 
been considered in Refs. [THUD]- It has been suggested 
that, since the conformal fluctuations formally give a neg- 
ative contribution to the Hamiltonian density, they can 
compensate an energy of the zero point fluctuations of 
matter fields, electromagnetic and gravitational waves. 
Nevertheless, as was shown in 14], it is not the case, 
at least for the gravitational waves. The reason is that 
the scale (conformal) factor is attributed naturally to the 
slowly varying background, and thus, it can hardly be 
modulated by gravitational waves. 

It should be added, that the decomposition of a metric 
into the conformal factor and a remaining part: — 
tt 2 g^v = e 2 "^!/ seems very natural as it was suggested 
by York [l5|, who discussed properties of the conformal 
geometry given by <? M „. 

In this paper we investigate the classical solutions of 
some simple models to determine a possibility of the scale 
factor modulation. This means that one explores some 
particular gauges and initial conditions, which result in 
an oscillating scale factor. This does not mean that the 
scale factor is an oscillating quantity under an arbitrary 



initial conditions. The aim of our work is to demon- 
strate, that the scale factor oscillations are not forbidden 
in general (e.g. owing to constraints). 

Although, we investigate local oscillations of the scale 
factor here, the models admitting global oscillations 
[l9j are also of general interest. 



II. STRING IN A BACKGROUND SPACE 

The first model considered is a string against a curved 
background. For this model we have built a solution for 
which the scale factor is connected rigidly with the mat- 
ter fields so that a mutual compensation of their fluc- 
tuations and, thereby, the solution of the cosmological 
constant problem look almost trivial. 

The Lagrangian for a string against a background 
space is given as [20| 
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where 4>(x, r) = {(/>i(x,t), ^(x, t) . . . ^n(x, t)} are the 
scalar fields, a(x, r) = In a(x,r) is the logarithm of a 
"scale factor" a(x), Af and Afi are the lapse and shift 
functions, respectively; M p is a constant (an analog of 
the Plank mass); primes denote the derivatives over r. 
This Lagrangian is not deduced directly from the gen- 
eral relativity action, however, it can be reduced to the 
form corresponding to a flat Friedman- Robertson- Walker 
(FRW) universe filled with the homogeneous scalar fields 
(i.e. if both a and <j) do not depend on x). Thus Eq. ([1]) 
can be considered as a simplest inhomogeneous general- 
ization of the FRW Lagrangian. 
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FIG. 1: Evolution of a'(x, t) at x — 1.5 according to Eq. (fTO) . 
Numerical values of parameters are y/N /M p = 1.5, A = 1.3, 
k — 2. Dashed line corresponds to < a' >, where the spatial 
averaging is implied. 

Varying over Af and Af\ yields the Hamiltonian and 
momentum constraints, which, in the case of Af = 1 and 
Afi = 0, take the form: 

-H = e 2a {-M 2 p a' 2 - M 2 p (8 x a) 2 + </>' 2 + (^0) 2 ) = 0, (2) 
V = e 2a (-M 2 a' x a + 0' d x 4>) = 0. (3) 

Also, one can obtain the equations of motion [20| 

<P" - d xx cf> + 2a'(p' - 2d x ad x <P = 0, (4) 
M 2 a" - M 2 d xx a + M 2 a /2 - M 2 (d x a) 2 

+4>' 2 - (d x <f>) 2 = o. (5) 

The constraints evolve according to the equations 

W = d x r, 

V = d x H. (6) 

That is, if the constraints equal to zero initially, they 
remain to be zero during all the evolution. 

Some particular solution can be easily found for the 
system @, ©, d3J, ([5]). Let us assume a(x,r) = 

$Mx>t) = ^fc(x,r) = ...^ n ( x ,t), then the 
equations of motion (0]),© are reduced to the single 
equation 

a" - d xx a + ^ (a' 2 - (d x af) = 0. (7) 

Also, one can be sure that the constraints ([2]), ([3]) are 
satisfied, i.e. T~L — 0, V = 0. 

It is interesting to note, that the analogous equation 
arises in the the Szekeres theory of colliding waves [2lj |. 
However, the initial conditions for this equation in the 
Szekeres theory are restricted by the constraints, while 
in the present theory the constraints are already satisfied 
and the initial conditions for Eq. ([7]) are quite arbitrary. 
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FIG. 2: Evolution of a'(x,t) (a), B'(x,t) (b) and V'{x,t) (c) 
at x — for the inhomogeneous Bianchi I model. Numerical 
values of the parameters are ao = —2, ny ~ 1/100, nso = 
2/100, Bo = 0, Vo(x) — cos x/ 100. Dashed lines correspond 
to the mean values < a' >, < B' > and < V' >. 



The oscillations against a background can appear if the 
initial values meet the condition 

|6( e a(a:,0)|»K(ar > 0)|. (8) 



In the opposite case, a monotonic evolution dominates. 
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Here, the global evolution is denned by the value of k 
and the oscillations are defined by the amplitude of A. 
The example of evolution at some fixed coordinate x is 
shown in Fig. Q] 

From the point of view of the "energy balance" spec- 
ified by the Hamiltonian constraint ((2]), the oscillations 
of scalar fields do not affect a slow "global" (i.e. "back- 
ground" ) evolution because they are compensated by the 
oscillations of scale factor a. 



ILL INHOMOGENEOUS BIANCHI I MODEL 



FIG. 3: Evolution of the dispersion V < a' 2 >. Dashed line 
corresponds to < a' >. 



Solution of Eq. ((7} can be written in the form of 

I \ M P 1 ( 1 #«o(^+r) . 1 ^K ao ( X ~r) 

alx, t) = — ==ln -e m p H — e M " 
V ' 2y/N \2 2 



~mZ 



x+t 



W , (9) 



As the next example, let analyse the inhomogeneous 
Bianchi I model. Hereinafter, we shall decompose the 
metric into the scale factor e 2a and the components of 
conformal geometry, and use the conformal time r. The 
argumentation for using the conformal time gauge is that 
it allows representing the equation of motion for a scale 
factor of universe in a form containing a difference of 
the potential and kinetic energies of the field oscillators 
0,11. 

Thus, the metric looks like 



where the functions ao(£) and fc(£) are connected with 
the initial values of a and its derivative so that a(x, 0) = 

ao(a:), a'(x,0) = e~ 2 ^ ao[x) k{x). 

As an example, we take k = const, ao(£) = ^4cos£. 
Then the solution of Eq. @ takes the form 



2 A cos 



a(x ' T)= ^ in (K 
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ds 2 = e 2a [ dr 2 - e~ AB dx 2 - e 2B+2 ^ v dy 2 



_ e 2B-2V3V dz 2 l (n) 



where the functions a, B, V depend on x and r. 

The Einstein equations allow obtaining the Hamilto- 
nian and momentum constraints as well as the equations 
of motion: 



H = -e 2a (-a' 2 + B' 2 + V' 2 ) 



1 



2 " li; ' \{d x a) 2 + ±d xx a + 7 -(d x B) 2 + ±d xx B + ~d x ad x B + \{d x V) 2 ) . ( 12) 



1 



1 



V = e 2a --d x aa' + d x BB' + -d x aB' + d x W + -d x B' + -d x ol 



n" < w { d xx a + {d x af + {d x Vf + 7 -{d x B] 2 + ^8 XX B + 4d x Bd x a^j + a' 2 + V' 2 + B' 2 = 0, 

B" + ie 4B (d xx B + 2(8 X B) 2 + 6(d x V) 2 - 2(d x a) 2 + 2d x ad x B + 2d xx a) + 2B'a' = 0, 
V" + 2W - e 4B {d xx V + 2d x Vd x a + Ad x Bd x V) = 0. 



(13) 

(14) 

(15) 
(16) 



r 



It should be noted that the equations of motion can be 
also obtained from the Hamiltonian H = J T~idx directly. 
The constraint evolution is given by © as it takes place 
for a string. 



Let ask: would the ripples of scale factor appear if it 
was initially homogeneous. Since there is no analytical 
solution, we shall investigate a particular solution of this 
model numerically. 
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The initial condition are taken in the form 

a{x, 0) = ao = const, 
B(x,0) = B = const, 
V(x,Q) = V {x) = Acos(x), 



V'(x, 0) = exp(— 2a )TT V (x), 
a'(x, 0) = exp(— 2cto)ir a (x), 
B'(x,0) = exp(-2a )ir B (x), 



(17) 



(18) 



where the velocities are expressed in the terms of initial 
momenta ttv(x), ir a (x), ttb(x), which are chosen as 

7Ty(a;) = const, 



°+ 4Q( %' 2 (x) 



2(7r BO - Wq{x)-k v ) 
ir B {x) = -w a (x) - 37ryVb(a;) + ir B o, 

ttbo = const, (19) 

i.e. in such a way that the constraints ([2]), (|3]) are satisfied 
exactly. 

The evolution of time derivatives at x = is shown 
in Fig. [3J One can see, that a'{x,t) evolves non- 
monotonically and can become negative at some mo- 
ment of time and in some spatial point x, i.e. the scale 
factor decreases here. Analogous oscillations appear for 
B' . Figure [3] shows the evolution of the dispersion of a', 
which finally becomes much greater than the mean value 
< a' >. 



theory. Let represent the variables as a sum of a spatially 
uniform part and a perturbation 

a(x, t) = a (r) + ai(r)e lkx , 
B(x,t) = B (t) + B x {r)e ikx , 
V(x,T)=V (T)+V 1 ( T )e ik *. 

(20) 

In the zero-order in perturbations, the equations of 
motion are 



V " + 2V 'a ' = 0, 
B " + 2B 'a ' = 0, 



0() 



" 1 Bn' 2 + Vn' 2 



a ' 2 = 0, 



(21) 




FIG. 4: Evolution of a'o (dashed curve), oscillations of a[ 
(solid curve), and V{ (dotted curve) according to Eqs. (12111 . 
121. 



IV. LINEARIZED THEORY FOR THE 
INHOMOGENEOUS BIANCHI I MODEL 

The origin of the oscillations demonstrated above can 
be easily explained within the framework of a linearized 



and the only constraint is 

U = \e 2a ° (B ' 2 + V ' 2 - a() ' 2 ). 
In the first-order in perturbations, one has 



(22) 



Vi" + k 2 e iB °Vi + 2V Q 'a 1 ' + 2V 1 'a ' = 0, 
Bi' + 2B 'ax + 2Bia ' - \k 2 e iB °{B 1 + 2ai) = 0, 

ax" + 2ao'a 1 ' + 2B 'B 1 ' + 2V 'V 1 ' + \k 2 e ABa {2B X + 3ai) = 0, 

Hi = -iB Q 'Bl' - 3aiS ' 2 + k 2 e iB °B 1 + k 2 e 4Ba ai - Wq'Vx - 3aiF ' 2 

-l-Sao'c^i' + 3aiQ!o' 2 = 0, 
V\ = —3ikBiB(/ — 2ikoi\B§ — ikB\ — HkViVo + ikaicto' — ikai' = 0. (23) 

I 



However, only one of the two last constraints turns out to be independent. Thus one may exclude B\ and after 
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FIG. 5: (a) Solid line is the mean value of the second-order contribution of £i to the Hamiltonian constraint. Dashed line is 
the background evolution i.e. Ho given by (|22[) . (b) Solid line is the mean value of the second-order contribution, but without 
taking into account a "pure" gravitational wave < H2V >• Dashed line is the contribution of a pure gravitational wave, namely 
< H 2 v >■ 



some computations comes to the equations for a%,Vi, which can be solved numerically: 
I 



ax 



Vi" + k 2 e iB °V 1 + 2V 'ai + 2ViW = 0, 
18B 'B "ai' + le SB °k 4 ai + e iBa k 2 (aVq (Vi - 3ViB ') - 4-Bo'ai' + a x {l (B ') 2 - 2B 

9 (Bn') 2 + e 4B «k 2 



r 



0. (24) 



Evolution of q;'i(t), V[(t) is shown in Fig. HJ which 
demonstrates that the amplitude of oscillations of a[ (r) 
increases with time. For the primary system (|12[) . (|13[) . 
(HU) , (ITS1) - (HU), the oscillations (Fig. 0} transform into 
the narrow spikes due to nonlinearity (see Fig. HJa)). 



V. RIPPLES FROM THE POINT OF VIEW OF 
THE ISAACSON'S THEORY 

The Isaacson's theory [24| considers an influence of the 
ripples on the evolution of background. An interpretation 
of the Isaacson's theory is quite straightforward and does 
not depend on details of the averaging procedure. Let 
expound its main aspects. Let suppose, that there is a 
set of equations 



A(e(x,r)) = o 



(25) 



for some quantities £, which have a strongly inhomoge- 
neous (oscillating) space-time behavior. Then, one may 
define an average quantity 



£o(t) =<£(t,x) >= - 



£(r, x)dx, 



(26) 



where we shall consider only spatial averaging for the 
sake of simplicity. 

As a result, the quantity £ can be separated as 



£(x,r) = £ (t) +£i(x,t), 



(27) 



where £i(x,t) — £(x,t) — £q(t). It is evident, that the 
mean value of ^i(x,r) equals to zero. 

Eqs. (|25|) can be expanded in powers of £i(x, t): 

A{t( Xi r)) = A (&(r)) + Ai(£ (t), ^ (x, r)) 

+A 2 (Ur)^i(x,r))+--- = 0. (28) 

The first statement of the theory considered reads that 
the spatially inhomogeneous part obeys the linear equa- 
tions A\ = 0, and, in fact, is equivalent to the lineariza- 
tion considered above. This statement is not valid in our 
case, because the ripples of B and a are strongly nonlin- 
ear. 

The second statement reads that the evolution of back- 
ground is determined by averaged terms of the second 
order in ripples: 

Ao(Ur)) = -< A 2 (Co(r),Ci(i,r)) > . (29) 
It is a consequence of the averaging of Eq. 
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Let consider only single equation: the Hamiltonian into the Hamiltonian constraint and obtain the second 
constraint (TT2"j) . We substitute order contribution in the form of 

a(x,r) = oo(t) + o>i(t,x), 
B(x,t) = B (t)+B 1 (t,x), 
V(x,t)=V (t) + V 1 (t,x) 

(30) 



U-2 



2a x d x B x B^ 
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{d x aiY 



-d x Bid x ai 



-aid xx a\ 



I (5, Si) 2 
6 

1 



-Bid xx ai 



+ \ (dA) 2 + 2a 1 d x V 1 V ' + ^ (V X 'Y - 2a 1 a 1 'a ' 



1 / a2 



■ ai 



'So') 



r 



(31) 



Using the results of the previous numerical calculations 
(Sec. one can see that the second statement of the 
Isaacson's theory (i.e. Hq as — < H2 >) is satisfied with 
a high accuracy as it is shown in Fig . [S] (a) . 

Let distinguish a "pure" gravitational wave contribu- 
tion 



1 



n 2V = e 2a °- mv) " + 



(32) 



in "H 2 - As one can see from Fig. [5] (b), the plot of 
< V.2 — 'H.2V > is, in fact, symmetric relatively the time 
axis. This means that if one performs the time averaging 
in addition, then one can conclude that < H2 — 'Hv2 > 
does not contribute "at the mean" to the evolution of 
background, i.e. the evolution of background is deter- 
mined by a "pure" gravitational wave < R 2 v >■ 



VI. THE YORK CURVATURE FOR THE 
NONUNIFORM BIANCHI MODEL 

Why is the V-variable considered as that correspond- 
ing to a "pure" gravitational wave? To answer this ques- 
tion, one may use the description of a conformal three- 
geometry in the terms of the York curvature [TEl [25j 

1 



Y 



ab 



(33) 



Other components of Y ab equal to zero. As can expect, 
the York curvature is conformally invariant, i.e. it does 
not depend on a. From Eq. (|34[) . one can see that V = 
results in Y = and the space becomes conformally 
flat. This justifies the consideration of the V- variable as 
representing a "pure" gravitational wave. 



VII. DISCUSSION AND CONCLUSIONS 

Let c omp are the model considered with the well-known 
Gowdy 26-29] model described by the metric: 

ds 2 = e- T - x ' 2 (e 4T dt 2 - dx 2 ) - e 2T (e 2 W + e^dz 2 ), 



(35) 



where T, A, /3 are the functions of t 1 x only. Denoting 

P = Vsv, 

T = B + a, 
A = 6(B-a), 



(36) 



we come to the metric 



ds 2 



„2a 



((e 2a dtf 



-4B 



dx 2 



_ e 2B+2V3V dy 2 
e 2B-2VSV dz 2) 



(37) 



where g ab denotes the spatial part of the metric (IT), Rab which I s exactl y the metric dl if the conformal time 



is the Ricci tensor of the three-metric g a i>, e ae f is the fully 
antisymmetric symbol and the ; e means the covariant 
derivative over the metric g a b- Direct calculation gives 
the following result 



dr 



s dt is used instead of t. 



Y 2 



= V3e 4S 



9d x Bd xx V + 3d x V[d xx B + 6(d x B 



2(d x V) 2 )+d xxx V . (34) 



For the Gowdy model, it is shown that the function 
T can be chosen as the function of the time t only (2(| 
H, and the function A ceases to depend on x at 
t — » 00. Thus, in principle, the ripples of a = \{T — 
■|) found above can be eliminated by some coordinate 
transformation. 

It should be noted, that, if the gauge invariance turns 
out to be violated in quantum theory, the ripples of the 
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scale factor in the conformal time can appear. Since, as 
it has been demonstrated above (Figs. [31 [5]), the ripples 
consist of the sharp spikes, their spectrum will be broad 
and "noise- like" . It should be taken into account, that 
such spikes have appeared when the Bianchi I model was 
considered. As has been demonstrated in Sec. IIV1 the 
unstable rising modes exist for this model. In contrast, 
it is known, that the isotropic FRW space-time is stable 
against small perturbations [30| . Thus, for the isotropic 
space-time, the ripples could exist only in a strongly non- 
linear regime and this issue needs a further analysis. 

Concerning a contribution of the ripples to the energy 
balance, one may conclude that such ripples occur in both 
string and Bianchi I models. But, in the string model the 
ripples of the scale factor and the oscillations of the mat- 
ter fields compensate mutually euch other and, thereby, 
the background evolution is not affected by them. For 
the Bianchi I model, the results of Sec. [V] demonstrate 



that the ripples of the scale factor compensate only the 
oscillations of the B-ficld and the background evolution 
"at the mean" becomes to be determined by the "pure" 
gra vitational wave V in agreement with the results of 
14|. Nevertheless, such a conclusion can become invalid 
in the quantum gravity case, since a strong gauge wave, 
consisting of the scale factor a and the £>-ficld could 
mixed up the gravitational waves and contribute to the 
background evolution, as well. 

It should be noted, that both models (i.e. the string 
and Bianchi I ones) have been considered on the classi- 
cal level that cannot shed a light on the solution of the 
vacuum energy problem. This results from the fact that 
a compensation of only zero-point fluctuations of matter 
fields is required in the quantum case. But in the clas- 
sics, it is impossible to distinguish the fluctuations in the 
ground state from those in an excited state. 
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